Fluctuations in quantum dot charging energy 
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We demonstrate that enhanced fluctuations in the charging energy of chaotic semiconductor 
quantum dots result from strongly scarred wave functions. We demonstrate that strong scars linger 
at the Fermi surface due to their enhanced Coulomb interaction between spin up and spin down 
states. We present density functional and spin density functional calculations to elucidate the 
temperature, dot shape and spin polarization dependencies of the fluctuations. 



A simple, physically appealing paradigm for the statis- 
tical behaviour of the ground state of chaotic, interacting 
systems does not yet exist. Based on the concept of level 
repulsion, the behaviour of eigenvalues and vectors of 
complex systems has been successfully treated for many 
years via random matrix theory (RMT) Q]. As recent 
experiments on quantum dots demonstrate however 
fluctuations in the ground state energy as a function of 
electron number N, are not well described by RMT since 
the required hypothesis of constant interaction (CI) ap- 
pears to be violated. On fundamental grounds this failure 
is to be expected since the ergodicity implicit in RMT is 
unconstrained by any minimization principle which must 
be an essential ingredient of any theory of the ground 
state. 

The charging energy Eq of a quantum dot, defined be- 
low in terms of the ground state energy at differing N, 
provides an important archetype of an interacting system 
insofar as particle number, disorder ||, the "chaoticity" 
j| and size can be varied. Attempts to understand the 
scale and behaviour of fluctuations in the charging en- 
ergy in quantum dots have included numerical studies 
of tight binding Hamiltonians for small particle num- 
ber as well as study of the statistical behaviour of 
electron-electron interactions in a circular, disordered dot 
based on the random phase approximation (RPA) for the 
screened interaction and a sigma model calculation for 
the eigenfunction correlator ||. 

Exact diagonalization studies such as in Refs. [ELpI give 
valuable insight into trends in the ground state energy, 
but they are necessarily limited to small particle number 
and in addition they depend on a reduction of a com- 
plicated set of matrix elements to a small set of energy 
parameters. Reference ||, on the other hand, considers 
explicitly direct and exchange Coulomb matrix elements, 
however the distribution of these elements are evaluated 
using a wave function correlator for a non-interacting dis- 
ordered system. Consequently there is no minimization 
principle and no self-consistency in the theory. Further- 
more, both the diagonalization approach and the method 
of Ref. H provide only statistical information on the 
charging energy. Neither method can realistically incor- 
porate a varying gate voltage and therefore the physical 
process whereby fluctuations occur is unapproachable. 



By contrast, a detailed study of the evolving, self- 
consistent eigenf unctions, free energy and level struc- 
ture of a realistic dot via density functional (DF) theory, 
which we present here, reveals the specific mechanism 
whereby Eq fluctuations deviate from the "RMT+CI" 
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FIG. 1. Ec vs. N and spectrum near Ef (lower two pan- 
els). Curves (a) and (b) differ from (c) by non-plunger gate 
voltages, hence varied pattern of peaks, (b) and (c) offset by 
0.04 and 0.08 Ry* resp. (a) and (c) are pure Hartree, (b) 
includes E xc though LDA. T = 0.1 K throughout. Spectrum 
corresponds to (a) and (b). Upper right panel shows levels 
from SDF calculation, same parameters as (b), arrows indi- 
cate Ec peaks. Note at N = 92 and 96 spin polarization of 
dot collapses. Inset: gate pattern and typical self-consistent 
effective 2D potential at 2DEG level. 

predictions. We show that large fluctuations in Ec are 
related to strong scars; remnants of periodic orbits in 
the classical confining potential. When located at the 
Fermi surface Ep these states, being quasi-one dimen- 
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sional, tend to create an inhomogeneous potential, in con- 
trast to other more "chaotic" states which occupy more 
smoothly the entire dot area. Consequently, when a scar 
is at the Fermi surface a gap is created to the next state 
which results in additional gate voltage spacing to the 
next Coulomb oscillation. Furthermore, we find that the 
direct Coulomb matrix element between a state and itself 
(up spin and down spin), denoted W pp = U p , is generally 
greater than between different spatial states W pq , p 7^ q, 
and that U p for scars is particularly large. Therefore, 
as we verify via full spin density functional (SDF) calcu- 
lations for certain parameter regions, the dot is almost 
always spin polarized 0. Further, scars tend to remain, 
half filled, at Ep (as V g is varied) due to their large U p , 
while other, more homogeneous states pass through the 
Fermi surface Q (in the process gradually screening the 
scar state) . Finally, when it becomes energetically impos- 
sible to draw down a still higher state to Ep, the filling 
of the second spin state of the scar does occur, leaving a 
gap to the next state, as mentioned above. At this point 
not only the scar, but all states are doubly occupied and 
spin polarization of the dot collapses. 

We calculate the self-consistent electronic structure 
Q for a small quantum dot with the wafer profile and 
gate pattern of the device used by Sivan et al. [Q. The 
nominal two dimensional electron gas (2DEG) density 
is 3.1 x 10 _ 1 1 cm" 2 . Typical values for the average level 
spacing (not assuming spin degeneracy) and charging en- 
ergy are A ~ 0.012 Ry* and E c ~ 0.18 Ry*, respectively 
(1 Ry* « 5.8 meV). The gate pattern is rectangular 
(inset Fig. 1) but due to excess metal in one gate the 
effective confining potential is classically chaotic, a fact 
which is substantiated by the statistics of the single par- 
ticle level spacings (not shown). As in the experiment 
we vary the "plunger" gate voltage V g (inset Fig. 1), 
typically in steps of 5 mV, the other gates being fixed 
to isolate the dot from source and drain. For each V g 
the electronic structure is calculated at three values of 
dot chemical potential Udot — 0.0, ±0.2 Ry* . States are 
filled (and N thereby determined) according to a Fermi 
function, which except for very small dots is expected to 
give a good approximation to the full grand canonical 
ensemble ||. Thus, in the calculation N is not nec- 
essarily an integer. The values of [idot are chosen to 
change N by about ±1, relative to the "equilibrium" 
value (/j,dot = pleads = 0). The free energy F(N,V g ) 
is calculated according to ref. j| and Ec is defined as 
d 2 F(N, V g )/dN 2 , which we compute discretely. 

Traces of Ec{N) (N from fidot = solution) are 
shown in Fig. 1. The quasi-periodic peaks are associ- 
ated with the onset (e.g. N — 87, 114) or completion 
(N = 92, 96, 102, 124) of scar filling, the former being 
accompanied by downward fluctuation of states below 
the scar state, the latter with gaps above the scar state. 
Changing "non-plunger" gates varies the precise pattern 
of oscillations, but not the basic structure. Varying a 



"back gate" (not shown) also produces qualitatively sim- 
ilar characteristics. 
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FIG. 2. Selected eigenfunctions (moduli squared) and grey 
scale of Coulomb matrix elements at N « 124 (cf. curves (a) 
and (b) Fig. 1). Inset gives close-up of spectrum from Fig. 1, 
scar state indicated by triangles remains near Ef- 

Standard DF theory is a spin degenerate mean field 
theory. Even when exchange-correlation E xc is included 
in the local density approximation (LDA) it is known that 
some uncorrected self- interaction remains ]l0| . The ten- 
dency (Fig. 1) for levels to cluster at Ep results from the 
energetic advantage of occupying several states partially 
rather than one or two states fully. At finite temperature 
T such partial occupancy can be justified as an approx- 
imation to the weighted ensemble of many-body states. 
Further, in this regime, Wigner crystallization is unim- 
portant. Nonetheless, a more precise description, where 
the role of spin polarization is made manifest, can be ob- 
tained via SDF calculations. As these are numerically 
more intensive we only calculate the behaviour for vari- 
ous representative parameter regions. Our main physical 
conclusion is that the electrostatically driven tendency 
for levels to cluster, partially filled, at Ep in ordinary 
DF theory is replaced by an electrostatically driven spin 
polarization in SDF (e.g. it is energetically favorable to 
occupy two different spatial states singly as opposed to 
occupying either of them doubly). 

The spectrum for the SDF case is illustrated in the 
upper right panel of Fig. 1. The Ec fluctuations (not 
shown) occur at the same N values as in the DF cal- 
culations and coincide with the arrows. Note that the 
fluctuation at N = 87, which corresponds to the filling 
of the first spin state of a scar, is qualitatively different 
from those at N — 92 and 96, where scars become doubly 
occupied. In the latter cases, since U p of the scar is much 
greater than that of other states, double filling generally 
only occurs when all other states are either doubly filled 
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or empty. Therefore |0, when this state fills, all states 
of the dot become identically spin degenerate. 

Figure 2 shows wave functions and Coulomb matrix 
elements in the vicinity of a typical Ec fluctuation, 
N sa 124. Coulomb matrix elements are calculated using 
the kernel of Poisson's equation from the self-consistent 
calculation. Thus screening by the gates is included au- 
tomatically HJl . We have identified four stable orbits for 
this device of which two are observable in Fig. 2. Here 
ri(32 = 2 (n p = occupancy of state p) and = 0. The 
grey scale plot of Coulomb matrix elements shows that 
this state, as well as p = 60, another scar along a slightly 
different orbit, have higher values of U p pj. Scars oc- 
cur in families as additional nodes are added along the 
same orbit (compare p = 62 and p = 50 and note their 
interaction). The Coulomb matrix elements given here 
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FIG. 3. Potential difference contours showing influence of 
p = 62, scar state filling (left) and that of more chaotic state 
p — 64 (cf. Fig. 2). Plunger gates indicated. 

are between self-consistent states fl4}| . Thus, even though 
Wq2, 63 is relatively small (~ 0.18 Ry*), it is the interac- 
tion of these two states which produces the gap at Ep 
and the fluctuation of Ec at this N. The strong 1-d scar 
p = 62 at Ep acts like a single atomic valence electron 
and strongly influences the potential (Fig. 3), creating 
a central "ridge" which modifies p — 63, reducing the 
Coulomb interaction between p — 62 and p = 63. 

A comparison between our calculations and the experi- 
mental results of Sivan et al. [^) shows a striking discrep- 
ancy. While the gate voltage spacing fluctuations (A2 
in Ref. (2)) are symmetric, the fluctuations of Ec com- 
puted here all proceed upward from a smoothly decreas- 
ing base. Furthermore, a comparison of the histograms 
of the fluctuations (not shown) reveals that the fluctua- 
tions in Ec have an r.m.s. deviation of only about 5% 
of Ec , a result that is smaller by a factor of two or three 
from the data for A2 . These discrepancies emphasize the 
importance, heretofore not fully appreciated, of distin- 
guishing between the "inverse compressibility" (what we 
call the charging energy) at fixed V g and the gate voltage 
spacing (normalized by the capacitance ratio a, see be- 
low) between Coulomb oscillations. These two quantities 
are equivalent only within the constant interaction model 
and, as we now show, small, upward variations in Ec can 
produce much larger, symmetric variations in AV^. 



To a reasonable approximation we can write 
F(N,V g ) = ±E c {Vg)N 2 + b{V g )N + c, i.e. F is quadratic 



in N at each V g . In the CI model, b(V g ) 



xV„ where a 



is the ratio of the dot-gate capacitance C c i g to the dot self- 
capacitance C, and Ec — jr + J2p e j> ( a * T = 0). C, a 
and e° (the single particle level energies) are all assumed 
to be independent of V g . AV g is determined by writing 
F(N, V g ) = F(N+1, V g ) and F(N+1, V' g ) = F{N+2, V g ) 
and subtracting, resulting in: 



eaAV a 



C 



£ N+2 



e N+l- 



(1) 



However when C, a and e° depend on V g , the procedure 
for deriving Eq. [I] leads to 



{N + l/2)e 2 

C{Vg) 



+e% +2 (V g i)-e% +1 (Vg). (2) 



When, in particular, C fluctuates appreciably from V g 
to VI the influence on AV g can be magnified, as illus- 
trated in Fig. 4. Here we have calculated F(N, V g ) on 
sufficiently fine grids to determine AVg numerically. 

To complete the picture we note that the frequency 
of scars in the spectrum, and hence the number of large 
oscillations in Ec, decrease with TV. This is somewhat 
evident in Fig. 1 in that the spacing of oscillations in- 
creases with N. The scars are remnants of closed orbits 
and are thus effectively one dimensional. Their density 
of states (DOS) decreases as 1/VE, in contrast to the 
constant DOS for the 2D dot as a whole. Thus the over- 
all frequency of scars in the spectrum will diminish. In 
addition, the extent to which a scar at Ep pushes up 
the ensuing state should also diminish with N. (Again, 
the oscillations in Fig. 1 for larger N appear to grow 
smaller). The reason is twofold. First, screening by the 
dot electrons increases with N. Second, the influence of 
the inhomogeneous potential fluctuation produced by the 
scar on the energy of a typical chaotic state should de- 
pend on the wavelength of that state A in comparison to 
the dot size, i.e. the kinetic energy cost of avoiding the 
scar decreases with A. 

When one considers larger dots with complicated scars 
occupying much of the dot area, the distinction between 
scars and non-scars, in terms of Coulomb matrix ele- 
ments, tends to disappear jl5|]. In this case the fluctua- 
tions in Ec will depend on the dot size, the fraction of the 
dot area occupied by the states near Ep, the screening 
length A sc and perhaps dot geometry. 

Sivan et al. found that while the fluctuations in eaAV g 
were much greater than the average single particle level 
spacing A, the temperature dependence of these fluctu- 
ations was, seemingly, governed by A/ksT H. We have 
argued that the fluctuations in Ec should decrease with 
N due to decreasing frequency of scars and their 



3 




i i r 

117 118 119 120 
electron number N 

FIG. 4. Charging energy and gate voltage spacing vs N. 
Non-plunger gates differ from preceding plots, giving large 
fluctuation at N = 118. While Ec fluctuation is purely up- 
ward, fluctuation in V g spacing is symmetric and substantially 
greater in magnitude. Inset: Ec fluctuation for various T. 
Fluctuations vanish far below fcsT = Ec- 

decreasing influence on neighboring states. We can- 
not say definitively, however, whether these fluctuations 
should scale with A (dot area) or with Eq (dot radius), 
though since A sc ~ A , the former seems more likely 
in the large dot limit. Empirically, however, the fluc- 
tuations in Eq are comparable to A and, as noted, the 
enhancement of fluctuations in AV^ (Eq. |2|) derives from 
V g dependence of capacitances and energy levels. Conse- 
quently, (Fig. 4 inset), the fluctuations essentially vanish 
by 600 mK (sw 0.75A) even though this is a factor of 
twenty smaller than Eq itself. 

Finally, in recent studies by Patel et at, citation in 0, 
no correlation was found to exist between peak heights 
and spacings. We propose that a stadium geometry, 
with quantum point contacts (QPCs) placed along the 
long axis, should indeed show this correlation. Certainly 
one stable periodic orbit will be that which simply oscil- 
lates along the central axis, so that these states will give 
greatly enhanced peak heights. Since, according to our 
discussion, scars associated with this orbit will float at 
Ep for certain ranges of V g , one should find correspond- 
ing groups of large Coulomb oscillations separated by a 
large V g spacing from neighboring oscillations. 

In conclusion, we have demonstrated that strongly 
scarred, quasi-ld states, by virtue of their strong 
Coulomb interactions, produce quasi-periodic fluctua- 
tions in the charging energy of chaotic quantum dots. We 
have shown that the scale of those fluctuations is consis- 
tent with the single particle level spacing, but that fluc- 
tuations in gate voltage spacings of Coulomb oscillations 
are, due to V g dependence of capacitances and energy lev- 
els, substantially greater and differently (symmetrically) 
distributed. Finally, we have proposed a simple exper- 



imental arrangement which could elucidate the role of 
self-consistency in chaotic systems discussed herein. 
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We note that the application of a magnetic field B, by re- 
ducing the effective dimensionality of all eigenstates also 
imparts structure to the Coulomb matrix elements. The 
influence of B on the spectrum can be seen in Ref. 3, 
Fig. 13. 

This is to be distinguished from the Coulomb matrix ele- 
ments, between bare levels, screened by the dot electrons 
themselves, as determined, for example, via RPA. 
For larger N one expects scars to be associated with 
groups of eigenstates rather than our small dot case 
where they are represented by single eigenfunctions. 
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